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Abstract. - Noninteracting fermions, placed in a system with a continuous density of states, may 
have zeros in the A'^-fermion canonical partition function on the positive real /? axis (or very close 
to it), even for a small number of particles. This results in a singular free energy, and instability 
in other thermal properties of the system. In the context of trapped fermions in a harmonic 
oscillator, these zeros are shown to be unphysical. Our results are also applicable to mean-field 
canonical calculations for fermions. By contrast, similar bosonic calculations with continuous 
density of states yield sensible results. 



Long back, Lee and Yang [Tl[2] pointed out that the 
onset of a phase transition could be deduced by studying 
the zeros of the grand partition function on the complex 
fugacity plane. As the number of particles goes to infin- 
ity in the thermodynamic limit, the complex zeros tend 
to pinch the real fugacity axis, signalling a phase transi- 
tion. Later, Fisher [3] found a similar behaviour of the 
canonical partition function Zjf(P) on the complex /3 (in- 
verse temperature) plane near a phase transition. Given 
a singlc-particlc partition function Zi(/3), Zjv(/3) may be 
calculated exactly for noninteracting bosons or fermions 
using recursion relations [4] . The single- particle partition 
function may arise from a one-body trapping potential, or 
may be the result of a calculation in a many-body prob- 
lem. In this context, Miilken et al. [5] have studied the 
pattern of Fisher zeros for trapped noninteracting bosons 
in a harmonic oscillator. As the number of bosons was 
increased, the real positive /3 axis tended to be pinched 
at T = Tc, the Bosc-Einstcin condensation temperature 
for spatial dimensions d > \. The behaviour of the heat 
capacity as a function of T also showed a sharp peak for 
large TV at T = T^- As expected, this was found both for 
the exact single-particle density of states of a harmonic 
oscillator (HO), and for the corresponding (asymptotic) 
continuous density of states. Miilken et al. [5] examined 
the distribution of zeros on the complex j3 plane to classify 
the order of the phase transition in finite systems. Simi- 
lar analyses were done with the Fisher zeros in interacting 
statistical models by Janke et al. [6l[7]. 



Noninteracting fermions trapped in a HO have not been 
studied in this context, presumably because no irregular- 
ity or phase transition is expected. However, a classical 
version of this formalism has been applied in a model for 
the thermodynamic properties of nuclear multifragmenta- 
tion in heavy-ion collisions [8]. In this paper, we are es- 
pecially interested in the effects of the Fermi statistics on 
the analytical properties of the canonical partition func- 
tion. In particular, we compare results obtained from the 
exact single-particle partition function from the discrete 
HO energy spectrum, with those obtained using a contin- 
uous density of states approximation. The latter is widely 
used in the grand canonical formalism. With the exact 
HO density of states, the fermionic system shows no ev- 
idence of any phase transition in any dimension, as ex- 
pected. With continuous density of states, however, we 
obtain some very surprising results, to be described be- 
low. We conclude that for fermions trapped in a HO, or, 
for that matter, confined in a box, these zeros are spuri- 
ous. Wc should mention here that in quantum field theory, 
lattice calculations of free Wilson fermions also manifest a 
real zero on the analogous inverse coupling axis if periodic 
boundary conditions are used [9]. This also gives rise to 
anomalous properties of the system [10] . Our paper may 
stimulate the search for spurious zeros of the canonical 
partition function in other areas of physics. They may 
appear when a continuous density of states is used that 
relaxes the restrictions in long-range correlations, or the 
use of improper boundary conditions. The long-range cor- 
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relations in our case arose from the Pauli principle, but in 
other cases may conceivably arise from long-range inter- 
actions. 

The iV-particle canonical partition function is generated 
from a one-body partition function Zi{/3) by using a re- 
cursion relation [3] 
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where a = I for bosons, and (—1) for fcrmions, and 
Zo{/3) = 1. The particles are (otherwise) taken to be 
noninteracting and identical, but correlations coming from 
quantum statistics are included. In this paper, for analyti- 
cal simplicity, the particles are taken to be noninteracting, 
and in a d-dimcnsional isotropic harmonic oscillator. The 
corresponding one-body partition function is given by 
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where we have not included the spin degeneracy for 
fermions. This when substituted in eq. ([1]), enables 

us to calculate Z^ili) exactly for noninteracting bosons or 
fermions occupying the discrete HO spectrum of energies. 
The single-particle density of states, g{e) may be obtained 
from a Laplace inversion of Zi(l3), and may be written 
as a sum of a smoothly varying part g(e), and an infinite 
sum of oscillatory terms [11]. For example, for c? = 2, 
^(^^(e) = I]^Q(n+ l)(5(e— (n + l)/icj), and may be written 
exactly as 



2^cos f 27rfc— ^ 
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To see how this comes about, note that 



Z\{fi) = -cosech 



^ OO 



^ (4) 

with no restriction on j3. On taking the Laplace inverse 
with respect to e, the single pole at /3 = for n = yields 
the first term on the right-hand side of eq. ^ , and the two 
poles for every integer n > at /3 = ±2imT/huj give rise to 
the oscillatory terms. The continuous part of the density 
of states above is e/(fia;)^, which we denote by 'g^'^\e). In 
thermodynamic calculations, it is the continuous density 
of states g^^) that is used widelj0. The smooth one-body 



^This smooth density of states could also be derived by a semi- 
classical expansion of Z\ {/3) in powers of h 1131 . Note that 

cosech(z) = \ z — ■ 

z 6 360 

Putting z = tkJl3/2, we get 



l^l < TT. 
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(Thomas-Fermi) partition function in d dimensions is 



(5) 



whose Laplace inverse yields the continuous density of 
states 

{d-l)\{hu:Y ■ 

The zeros of Zn{I3) on the complex /? plane may be cal- 
culated for bosons as well as fermions. Note that the ana- 
lytical expression for Zm (/3) is determined by the choice of 
Zi{l3) given by eq. ([2]), or its continuous counterpart Zi{f3) 
defined above. Schmidt and Schnack [M] have expressed 
Zn(P) in terms of polynomials in exp(— fia;/3), which make 
the calculation of the zeros for the discrete cases manage- 
able. As stated earlier, the bosonic calculations had been 
done in Ref. [S] , and we repeated some of these to check our 
programs. Our calculations showed a clear tendency, even 
for as few as = 10 bosons, for Fisher zeros to approach 
the real positive /3 axis for d ~ 2, and more so for d = 3. 
This is shown in fig. [1] Although the actual distribution 



Im(|3) 



Bosons 
d=3 



- 

X if X 
X ii « 

xV* ' 
x^x"* 
" xx Vx 



A'= 5 
^=10 



-0.6 



-0.4 



-0.2 



0.2 
Re(P) 



0.4 



0.6 



0.8 



Fig. 1: (Color online) The zeros of the exact Zn{P), for N — 5 
and N — 10, obtained from Zi{l3) of eq. ((5)) as a function 
of /3 = 1/kT for noninteracting bosons in a three-dimensional 
harmonic oscillator well. 

of zeros on the complex plane differed for the discrete and 
the continuous density of states, both showed the expected 
behaviour as N was increased. These bosonic results have 
been studied [5], and will not be displayed further. 

We now present the fermionic results. Consider first 
the discrete density of states, which gives rise to Zi{[3) 
of eq. ([2]), and the resulting Zn for fermions. In fig. [21 
its zeros on the complex plane are displayed for N = 5 
and = 10 fermions in two and three spatial dimensions. 
There are no zeros other than the origin in one dimension. 
Since Z^IP) is known, the thermal properties (for real /3) 
are easily calculated. These show the expected regular 
behaviour. 



with the restriction that \fkij0\ < 2tt for the asymptotic series to be 
valid. On Laplace inverting with respect to /3, the smooth is 
again reproduced, followed by a delta function and its derivatives at 
the origin. 
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(a) Zeros of the partition function in two dimensions. 
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(b) Zeros of the partition function in three dimensions. 

Fig. 2: (Color online) The zeros of the exact Zn{I3), for N — 5 
and A'^ = 10, obtained from Zi{l3) of eq. Q as a function of 
P = 1/kT for noninteracting fermions in a harmonic oscillator 
weU: (a) d = 2 and (b) d = 3. 



When, however, Zm{I3) is calculated using Zi{l3) given 
by eq. ([5]) for the continuous density of states, we get very 
unexpected results. For even N , we find that one zero falls 
precisely on the positive real j3 axis, even for very small 
particle number (i.e., N = 2), regardless of the dimension 
d . We investigate the thermal properties of such a sys- 
tem when calculated using the canonical formalism. The 
resulting heat capacity is negative over a range of tem- 
perature, indicating instability of the system. This will 
be discussed shortly. In fig. [3l we display the zeros of 
Zn{P) on the complex /? plane for = 10 in one, two 
and three spatial dimensions. The associated heat capac- 
ities for these cases are also shown. Similar behaviour is 
also found for very small to very large even N. For odd 
A'', however, Zm{P) from the continuous density of states 
has no real zero on the positive /3 axis. This could be 
easily checked analytically, for example, for A^ = 3. Nev- 
ertheless, the complex zeros of Zn{13) get dense near the 
positive /3 axis and give anomalous thermal properties. 
This is displayed in fig. |4] for A^ = 15. Similar results are 
obtained for smaller as well as larger odd A^. 

The present calculations with the continuous density of 
states give unphysical results for fermions in a harmonic 
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(a) Zeros of Zjv (/9) for 10 fermions calculated using Zi {j3) given 
by cq.lO in three different dimensions. 
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(b) Heat capacity as a function of temperature for a 10- 
fcrmion system, under the same approximation as in (a) 
above. 

Fig. 3: (Color online) The zeros of the canonical partition func- 
tion as a function oi (3 = 1/kT (a) and the heat capacity as 
a function of t = 1/(3 (b) for a system of 10 noninteracting 
fermions with a continuous density of states, given by eq. ([6]). 



potential. For example, take the case of A^ = 2. From 
Eq. we get 



Z2(/3) = (Z2(/?)-Zi(2/3))/2 



(7) 



For the continuous density of states, replacing Zi{/3) by 
Zi{/3), and using Eq. dSJ for d = 2 with {huj = 1), we 
obtain 

Z2(/3) = (I//34 - l/(4/32))/2 , (8) 

which has a real positive root for /3 ~ 2. This behaviour for 
trapped particles is spurious, and can be seen as follows. 
Consider the single-particle canonical partition function of 
a one-dimensional system. The d-dimensional counterpart 
may be obtained by taking its d^^ power. Writing Zi{/3) = 
exp(-/3ei), we find, using Eq.Q 
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(a) Zeros of the partition function of 15 fermions for three dif- 
ferent dimensions. 



1000 



800 



600 



400 



200 




(b) Heat eapaeity of a 15-fermion system. 

Fig. 4: (Color online) The zeros of the canonical partition func- 
tion as a function of /? = l/ksT (a) and the heat capacity as 
a function of r = 1//3 (b) for a system of 15 noninteracting 
fermions with a continuous density of states, given by eq. ((6]). 



The second term on the RHS of Eq. ([7]) subtracts out the 
contributions where the two particles occupy the same or- 
bital. The resulting expression for Z2 (/?) in Eq. © above 
is in accordance with the Pauli principle, and cannot be 
zero for real positive /3j^ 0, if all states have discrete eigen- 
values. The roots of ^2(/3) = found from Eq. ([5]) may 
not be physically meaningful. This analysis may be ex- 
tended for larger (finite) N. Figure [S] shows the variation 
of entropy S with temperature for = 2, = 10 and 
A^ = 50 fermions trapped in a d = 2 harmonic oscillator, 
using the continuous density of states given by Eq. ([B]). 
The behaviour in all three cases is the same, and is dom- 
inated by an unphysical pole in S due to the anomalous 
Fisher zero. Note that the anomalous zero shifts from 
T = 1/2 to about T = A as N varies from 2 to 50. For 
N ~ 2, the anomalous zero is at T = 1/2, and near this T, 
(4r2_i/2)/(T2-l/4). This causes S'(T) to plunge to 
large negative values for T < 1/2, giving rise to a discon- 
tinuity at T = 1/2. This explains why the heat capacity 
is negative in this range of T. 

The energy fluctuation in a canonical ensemble is 
{{E^) - {Ef) = T'^C^T), where (...) denotes an en- 
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Fig. 5: (Color online) Entropy S for 2-dimensional fermion sys- 
tems with continuous energy-density spectrum. N = 2 (anal.) 
refers to entropy obtained from Eq. ([8]). 



semble average. Since the energy fluctuation may also be 
written as {{E— (_B))^), we normally expect this to be pos- 
itive. The change in sign can, however, be explained by 
examining the inverse Laplace transform of Eq. (|5]), which 
gives the two-fermion continuous density of states. This 
is not positive definite, and ensemble averaging (...) is no 
longer guaranteed to yield a positive definite answer. This 
is easily checked, for example, by calculating, for E > 0, 
the ensemble average (E) as a function of T. 

Quite generally, the exact density of states for particles 
moving in a potential, consists of a smooth continuous 
part, and a sum of oscillating terms [H]. Equation ([3]) 
for the HO is a specially simple example of this. It is 
the oscillatory contribution that gives rise to the discrete 
delta-function spikes to the density of states. The oscil- 
latory part has a large effect on the analytical behaviour 
of Zn{[3), specially for fermions. The fermionic Z^^P) 
is more sensitive since the recurrence relation given by 
Eq. ([T|) has alternating positive and negative terms. This 
is unlike the situation in the grand canonical partition 
function, given by 



lnZG(a,/3) 



g(e) In (1 + exp(a - /3e)) de , (10) 



where a — /3fi, and fi the chemical potential. It can be 
shown that the contribution of the oscillatory terms is 
exponentially small with a 1//3 dependence in the expo- 
nent |15] . The grand canonical formalism generally gives 
sensible results using the smooth single-particle density of 
states g{e). We show that the same 5(e), when used recur- 
sively to generate the canonical Zpf{l3), yields anomalous 
results subject to misinterpretation. 
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